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Scattering  Properties  of  a  Moving  Dielectric  Object 


FOREWORD 


This  report  deals  with  the  scattering  properties  of  moving  dielectric 
non-spherical  objects,  a  subject  that  has  been  treated  in  strict  collaboration 
with  profs.  F.  Borghese,  P.  Denti  and  G.  Toscano,  University  of  Messina,  Istitu- 
to  di  Struttura  della  Materia,  and  with  dr  O.I.Sindoni,  CRDC,  Aberdeen  P.G. ,  Md. 
USA. 

The  research  reported  in  this  document  started  in  November  1983  and 
is,  at  present,  under  further  development.  In  fact,  the  subject  is  rather  promi¬ 
sing  anc  far  from  exhausted  so  that  I  plan  to  continue  along  the  lines  sketched 
below  even  after  expiration  of  the  present  contract. 

In  view  of  the  vastness  of  the  subject,  I  will  presently  deal  with 
two  particular  aspects  of  the  scattering  from  moving  objects.  First  of  all  I  will 
consider  the  effect  of  the  change  of  shape  of  the  scatterers  as  a  consequence  of 
mutual  collisions.  This  model  is  that  commonly  used  to  describe  chemical  reactions 
and  is  thus  of  wide  interest.  The  second  aspect  I  will  deal  with  is  that  of  the 
scattering  properties  of  particles  of  very  high  and  very  low  refractive  index, 
with  possible  onset  of  resonance  scattering.  As  will  be  seen  below  this  phenome¬ 
non  has  been  widely  studied  in  the  past  more  as  a  mathematical  curiosity  than  as 
a  true  physical  effect.  However,  as  a  consequence  of  general  covariance  of  the 
Maxwell  equations  under  holonomics  space-time  transformations,  it  may  happen  that, 
a  fast-noving,  scattering  particle  may  appear  to  an  observer  in  the  laboratory  as 
an  object  of  very  high  or  very  low  refractive  index.  Although  the  investigation 
of  this  particular  subject  is  far  from  complete,  I  will  present  here  some  prelimi¬ 
nary  results  which  may  be  interesting  for  other  researchers  in  the  field. 


PART  A 


| 

EFFECT  OF  THE  CHANGE  OF  SHAPE  OR  OF  RECOMBINATION  OF  SCATTERERS  ON  THE  OPTICAL 
PROPERTIES  OF  A  POLYDISPERSE  MODEL  AEROSOL 

Let  us  consider  a  homogeneous  medium  within  which  a  number  of  scatte- 
rers  are  uniformly  dispersed.  The  scatterers  need  not  be  ail  identical  to  each 
other  but  may  be  of  several  kinds,  as  is  the  case  in  polydisperse  aerosol.  I  also 
assume  that  the  scatterers  are  nonspherical  but  that  their  individual  scattering 
properties  can  be  well  approximated  by  modelling  them  as  clusters  of  spheres.  The 
collective  optical  properties  of  the  dispersion  can  be  calculated  by  studying  the 
features  of  the  coherent  propagation  of  the  electromagnetic  field.  Strictly  spea¬ 
king,  to  be  completely  correct  one  should  include  in  the  equations  for  the  propa¬ 
gation  the  effects  of  multiple  scatterings  among  different  scatterers  and  the  ef¬ 
fect  of  motion  of  the  single  scatterers.  However  the  former  effects  can  be  negle¬ 
cted  if  one  restricts  itself  to  dispersions  of  low  density  as  is  the  case  of  actu 
al  aerosols;  as  for  the  second  effect  it  is  well  known  to  produce  broadening  of 
the  peaks  in  the  spectra  of  the  single  scatterers  and  no  further  detectable  effects 
in  view  of  the  low  speeds  attained  by  the  scatterers  in  actual  aerosols.  Anyway 
the  effect  on  the  collective  spectrum  of  the  dispersion  is  further  smooted  on  acco 
unt  of  the  random  and  uncorrelated  motion  of  the  scatterers. 

As  a  consequence  of  the  above  consideration  our  effort  has  been  directed  to  see  if 
some  effect  is  detectable  when,  as  a  consequence  of  their  motion,  two  or  more  scat 
terers  collide  and  recombine  to  give  scatterers  of  different  structure.  As  will  be 
seen  in  the  following  paper,  which  is  the  final  version  of  a  paper  enclosed  as  a 
part  of  the  Second  Periodic  Report,  and  that  has  been  accepted  for  publication  in 
II  Nuovo  Cimento  B,  the  effort  has  been  successful  for  the  change  of  shape  is  quite 
detectable  in  spite  of  the  simplifications  we  assumed  to  reduce  the  computations. 

Before  closing  this  section  I  want  to  stress  that  the  content  of  the 
above  mentioned  paper  has  been  presented  at  the  1st  Int.  Conf.  on  Aerosol  Res.  and 


Tech,  held  in  Minneapolis,  Mn.  USA  from  17  to  21  September  1984. 

Furthermore  I  attended  to  the  4th  Annual  Meeting  of  the  sector  "Spectroscopy" 
of  the  Gruppo  Nazionale  di  Struttura  della  Materia,  held  in  Alghero,  Sardinia, 
from  11  to  13  October  1984.  In  this  meeting  it  emerged  a  new  interest  in  both 
the  experimental  and  theoretical  aspects  of  physics  and  chemistry  of  clusters 
by  the  comunity  of  italian  researchers.  This  circumstance  and  the  results  ob¬ 
tained  in  the  past  incourages  me  to  continue  the  research  outlined  above  even 
after  the  expiration  of  the  present  contract. 
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Summary.  —  Id  a  preceding  pajter  tb«-  absorption  coefficient  foi  tbe 
propagation  of  tbe  electromagnetic  field  through  a  lot*  -density  ditj^-r- 
sum  of  nonspbencal  scattcrer*  of  random  oriental  ton  vat-  calculated 
by  modelling  tbetn  as  clusters  of  spheres.  In  tins  ]..i}>ei  tbe  spbcief  are 
allotted  to  rearrange  among  tbe  clusters  to  simulate  cbcmir.a]  reactions 
and  tbe  changes  inducted  on  tbe  absorption  coefhcienl  of  tbe  dis]>ersiOD 
are  ..discussed  in  ~a  few  significant  case*. 

PACS.  42.20.  -  Projiagatiou  and  transmission  in  inhomogeneous  media. 


1.  —  Introduction. 

It  is  well  established  through  experimental  measurements  that  the  refrac¬ 
tive  index  and  the  absorption  coefficient  of  a  dispersion  of  identical  scattered 
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M-nsibly  depend  on  their  shape  (*  •).  In  order  to  eompare  the  exj><  rimental 
leMillf  (or  tW  piopugilion  through  #ct»*l  imdiit  %  it h  the  theoretic.*!  predic¬ 
tions,  1  v\ <>  uppro.ic.he.*  are  mainly  used:  the  effective  dielectric-constant  theo¬ 
ries^*)  (t.p.  t In*  Maxwell-Garnett  theory  and  the  effective  medium  theory) 
and  the  coherent  field  theories  ("‘  (t.r.  the  calculation  of  tlx-  coherent  prop¬ 

agation  of  the  electromagnetic  field  through  an  a-*embly  of  scatterers).  Even 
the  effects  of  the  size  distribution  (:-’),  of  aggregation  of  particles  (*- **)  and  of 
the  orientation  of  the  scatterers  ( 11  M)  were  studied  through  these  approaches. 

For  what  concerns  the  theoretical  investigations  several  author.*  dealt 
in  different  ways  with  the  scattering  by  individual  nousphorical  objects  (“  *•). 
In  particular,  the  pie**-nl  authors  (’■)  a*  well  a*  Gi:ka1;i>t  and  Ai'sloos  (1(  u) 
proposed  independently  to  model  some  kind  of  nonspherica!  scatterers  as 
clusters  of  spheres  The  individual  optical  rt  - pease  of  these  model  tcaHerers 
etui  be  effectively  calculated  as  a  function  of  their  orientation  with  respect  to 
the  incident  field.  Furthermore,  in  the  framework  of  the  coherent  field  theory 
we  recently  devised  { 1  * )  an  original  formalism  intended  to  study  the  effects 
of  the  any  sot  ropy  of  the  statU-Tcrs  on  the  propagation  of  the  electromagnetic 
field  through  a  low-density  random  distribution  of  clusters  of  spheres.  In  our 
formalism,  the  choice  to  use  clusters  of  spheres  as  model  stutterers  allowed  us, 
by  the  use  of  their  specific  features,  to  study  the  propagation  through  an  as¬ 
sembly  of  stutterers  having  a  high  degree  of  complexity.  On  the  other  hand, 
the  randomness  of  the  distribution  and  the  choice  of  the  low -density  limit, 
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letting  u-  tli -Tt^.tT  il  belli  the  stntMical  averaging  pro<  e.luie-  and  j},t-  elb-cl 
of  uiulli]ilc  m  a*  icring  among  dilb-ienl  cliKicr.s,  1  m •n<-l  1 allowed  us  to  evi¬ 
dence  the  an v. sol  ropy  effects  in  which  wo  arc  interested.  In  .sj.il*-  of  the  j.re- 
vuh-nl  interest  in  metallic  Mintlcrcrs  {"’ :s),  we  revived  to  include  inlo  the  c)u-.iers 
only  dielectric,  i.c.  nonconducting,  kj> lures  with  frequency-independent  re¬ 
fractive  index;  in  this  way  the  efTect'  of  the  nonsj.heric.al  shape  of  the  individual 
dusters  aro  not  masked  at  all  and  prove  to  be  quite  important. 

Since  our  formulism  of  ref.  (,;)  also  applies  to  disje-rsions  containin':  more 
than  one  kind  of  identical  scalterors  [c.g.  polydisporse  aerosols),  the  duster 
model  is  very  suitable  to  describe  the  effects  induced  on  the  scattered  field  by 
modifications  of  the  individual  scatierers.  Such  modifications  can.  in  turn, 
be  described  in  the  framework  of  the  model  by  rearrangements  of  the  spheres 
either  within  the  same  cluster  or  among  different  clusters.  This  mechanism 
could  be  identified  with  thal  of  the  chemical  reactions  as  long  as  a  cluster  of 
spheres  cun  represent  a  molecule  front  the  point  of  view  of  classical  optics. 

In  this  paper  we  intend  to  show  the  effects  of  the  above-mentioned  «  chemical 
reactions  »  on  the  absorption  coefficient  of  the  dispersion  of  scatierers.  There¬ 
fore,  in  the  following  section,  we  concisely  report  our  formalism  for  the  ab¬ 
sorption  coefficient  of  a  low-density  random  dispersion  of  clusters.  In  the  last 
section  wo  apply  the  formalism  to  the  dispersions  of  the  clusters  modelling  the 
*  reagents »  and  the  *  products  »  of  a  few  sample  «  reactions  »;  then  we  collect 
the  results  of  the  calculations  so  that  the  effects  of  the  «  re.»otions  »  are  prof¬ 
itably  displafed.  ft 

We  conclude  these  introductory  remarks  by  observing  that,  even  in  this 
paper,  we  consider  only  clusters  of  dielectric  spheres,  i.c.  nonconducting,  with 
frequency-independent  refractive  indices  because  we  want  to  single  out  the 
influence  of  the  structural  modifications  of  the  scatterers:  in  any  case  our 
formalism  can  be  applied  to  clusters  of  conducting  spheres  witL  refractive 
indices  obeying  some  dispersion  relation. 

2.  -  The  constants  of  coherent  propagation  through  a  dispersion  of  clusters. 

referring  el>ewhero  for  further  details  (1:),  here  we  i\>:.t.iM]y  icc.il!  how 

the  refractive  index  and  the  absorption  coefficient  of  a  dispersion  of  identical 
clusters  within  a  matrix  can  be  calculated.  In  the  formalism  reported  below 
the  dispersion  of  clusters  is  assumed  to  be  random  and  its  number  density  low 


(*•)  J.  I).  En.KSOLr  and  II.  I*.  liRom*:  P/*y».  Bet.  B,  15,  lt44  11977). 

(*')  C.G.  Gk*m>vvi5i,  N.  Can  s  ulk  and  0.  IlrsDERl;  Solid  Stair  Commun.,  31,  248 
(1979). 

(”)  F.  CLARO:  Phys.  Err.  B,  25,  7875  (1982). 

(u)  E.  Kims:  Phyt.  Bee.  B,  26,  3440  (1982). 
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enough  to  allow  us  lo  di-icg. ird  tlio  mull iph-scniu-ring  ofloot-  among  dif- 
fcn-nt  dusters.  We  also  assume  that  the  sphere*  in  the  clusters  as  well  as  the 
matrix  includin':  them  an-  made  of  isotropic,  homogeneous  dielectric  materials. 

Let  us  consider  wlialever  low -density  dispersion  of  scatterer*  and  send 
through  it  the  circularly  polarized  plane- wave  field: 


Er  —  \/-  Eu  ur.  ex p  \ik  ■  r]  ,  \Br  —  >;£* , 


w  here 


ur.  —  (t,  -j-  irit,)l\/2 

with 

C„/£„  =  A„„ 

and  i/  —  i  1  according  to  the  helicity  (“).  The  (complex)  refractive  index 
matrix  of  the  dispersion  is  consequently  given  by  (IS) 

whore  the  index  »•  number?  the  soatterors  in  the  volume  T  and  f  is  related 
to  the  normalized  for  ward-scattering  amplitude  of  the  >-th  scattorer  through 

/  =  «?*■  f  ■ 

The  absorption  coefficient  and  the  customary  refractive  index  of  the  dispersion 
are  given  by  (”) 

—  21:  Im  (~A  rr ) ,  fh:  —  Ke  , 

respectively. 

Now  jre  recall  that  ,  if  independent  of  the  location  of  the  v-tb  scatterer, 
but  does  depend  on  its  orientation  with  respect  to  the  incident  held.  There¬ 
fore,  if  the  low  density  of  the  dispersion  is  not  too  low,  we  can  substitute  /,  ,(0) 
for  /  .  and  write 


—  Am  ■+•  (0)/m(6)dfi  » 


where  the  argument  6  is  a  collective  symbol  for  the  Euler  angles  that  individuate 
the  orientation  of  the  scatterers  and  A(0)  is  the  number  density  of  the  scat- 
terers  with  orientation  b.  At  this  stage  we  specify  the  soatterers  in  the  dispersion 


f**)  J.  D.  Jacksok:  Clot  tiro!  Elretrodynamir*  (New  York.  N.Y.,  1975). 

(“)  K.G.  Newton :  ScaUrring  Tttrory  of  IToirf  and  FariirU*  (New  York.  N.1  ./  1966). 
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as  clu-<t<Tc  of  spheres  and  assume  that  t h**i r  oriental ion  is  random.  Then  we 
express  j,  ,(0)  in  terms  of  the  amplitude.-  of  tin-  mull ijmlt-  field-  scattered  by 
any  duster  with  orientation  0  and,  nukin':  full  use  of  tin-  features  of  the  cluster 
model,  wo  Kuco«vdod  iu  factorizing  this  quantity  into  a  part  depending  only 
od  t lit*  structure  and  a  part  dependin':  only  on  the  orientation.  This  allows 
us  to  perform  the  integration  over  the  orientation  analytically  with  the  result 


(1)  -^ni;  —  «  ~ I-  4^"»j  ^  tu  X  ~r  »/  )  • 


For  later  purposes  we  rewrite 

( -  )  U'l  —  <Vl!  4"  X  1/  1 

i 

where  we  define 

^  =  +  l  I  TF,''1J'Tr’,t1'  ’ 

*W  =  f,X  4-  . 

K  *• 

t 

In  the  above  equations  3’  is  the  Dumber  density  of  the  cdustors;  the  quantities 


where  the  X1Jf’s  are  vector  spherical  harmonics  (*•),  are.  proportional  to  the 
multipolar  amplitudes  of  the  incident  plane-wave  field;  the  matrices  V'A  aDd 
U'*',  whose  meaning  will  be  discussed  below,  Lave  a  rather  complicated  ex¬ 
pression  for  which  we  explicitly  rely  upoD  ref.  (**)• 

By  their  very  definition  the  IT’s  and  consequently  the  depend  only 
on  the  direction  of  the  wave  vector  of  the  incoming  field.  The  P’s  as  well  as 
the  S’s  depend  od  the  magnitude  of  the  wave  vector  of  the  plane  wave  and  on 
the  refractive  indices  of  the  spheres  constituting  the  (identical)  clusters  in  the 
disjtersion.  The  l’’s  account  also  for  the  structure  of  the  clusters  and  need 
to  be  calculated  once  for  all  referring  to  a  cluster  whatever  with  arbitrary 
orientation.  As  a  matter  of  fact,  the  V' s  are  responsible  for  almost  all  of  the 
amount  of  computations  required  by  eq.  (1)  because  their  calculation  implies 
the  inversion  of  a  matrix  of  order  23 ',LU  {L„  -f-  2),  where  N,  is  the  number 
of  spheres  in  a  cluster  aDd  L„  is  the  least  value  to  which  the  sum  over  L  must 
be  extended  to  obtaiD  well-converged  results.  The  convergence  of  the  multipolar 
expansion  on  which  eq.  (1)  is  grounded  Las  beeD  discussed  elsewhere 
of  course  we  shall  return  to  this  subject  wheD  considering  the  specific  com¬ 
putation  reported  iD  the  Dext  section.  Finally  it  is  useful!  to  recall  that,  if  the 
clusters  possess  symmetry  elements,  such  as  reflection  planes  aDd  symmetry 


I  O  1 


j...\ : 


Iii'-'  avo,  I  V.  lit  v  ) 


axi-r,,  group  theory  cun  substantially  help  to  compute  the  U't  for  it  implies 
the  inversion  of  matrices  of  smaller  order  than  otherwise 


3.  -  Changes  in  tlir  absorption  coefficient  caused  bv  «  clieinieal  reactions  ». 

In  tho  case  of  a  low-density  dispersion  containing  more  than  one  kind  of 
clusters,  eq.  (li)  becomes  (,s) 


< ,  = 


\yf, 

a ,  l 


whore  the  index  a  refers  to  the  a-th  kind  of  clusters  with  number  density 
A  ^  - —  A  Co  • 

Now,  it  is  well  understood  that,  if  a  •  chemical  reaction*  rearranges  the 
spheres  among  the  clusters,  it  determines,  starting  with  the  relative  concen¬ 
trations  Co  of  the  *  reagents  »,  the  relative  concentrations  df  of  the  «  reaction 
products »  as  well  as  those  of  the  «  residues  of  the  reaction  »,  c'  (at  least  one 
of  these  latter  is  zero  when  the  «  reaction  *  is  over).  Therefore,  we  can  cal¬ 
culate  the  absorption  coefficients  before  and  after  the  *  reaction  *  and  denote 
thorn  by 


y'"=  N  Y  c 

/  —  a 


and 


respectively,  where 


4'=*(2'X,+  2VJ. 

*  P 


Clearly,  a  very  useful  quantity  in  detecting  the  Etructural  effects  on  the  ab¬ 
sorption^ 


_PY?- 


(1 


=  pA„  -f  1  , 


where  p  is  tbe  completion  index  of  the  « reaction  »  (0<p<l),  and 


A, 


Yl 

rf 


i. 


In  the  following  part  of  this  section  we  present  the  results  for  a  few  «  chemical 
reactions  ».  Before  any  explanation  and  comment,  however,  we  want  to  stress 
a  point  to  avoid  any  substantial  misunderstanding  of  tbe  purposes  of  the 
present  wort.  Our  results  refer  to  sample  calculations,  indeed,  and  do  not 


1  \  I  i  rj  (•)  *3  li}  •  il.i 
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pretend  to  give  a  realist  it-  ih-t  riptiun  of  aelunl  optical  plicnorm-na  :  Ion  many 
details  should  he  added  to  the  \ cry  st liem.il i»-  olust»  r  model.  Con-iqueutly 
ooe  should  not  ht.lher  neither  about  the  choice  of  tie-  ]»:irauie.ter«  defining  the 
structure  of  the  clusters  i.or  about  the  possibility  of  the  ♦  reactions  »  in  the 
assumed  low -density  limit  of  the  dispersion.  We  confirm  that  our  question  is: 
if  the  structure  of  the  cluster  changes,  how  much  is  the  absorption  coefficient 
am'clcu*  "1  in-  results  we  aie  going  to  diMUss  111. the  u-  answel  that.  1L  ge.liei.il, 

these  changes  of  structure  modify  to  a  significant  extent  the  coherent  prop- 
agation  of  the  field  through  the  dispersion. 

The  *  reactions  »  for  which  f,,  was  actually  comjiutod  are  listed  in  table  1. 


TaKLI:  1  -  Lift  ot  tin  t  rtartion*  »  tonnidiuJ  in  thin  y-uprr.  The  refractive  indices 
and  radii  of  the  elements  A.  It  and  C  are  tiA  =■  1.30,  «„  =  1.50.  r.r  =  1  40  and 
*,  =  l.Ou,,  rh  =  O.SOu,,  rc  =  o.Tati,.  where  ti,  is  an  arbitrary  unit. 


K tact  ion 

Example 

1) 

3B,  —  2r’, 

30,-20, 

2) 

2AB  4  B,  —  2AB, 

2C04  0,  —  2C0, 

3) 

2AB,4B,—  2AB, 

2S0,4  0,  —  2S0, 

4) 

A, -33,  —  2AB, 

N,4  3H,  — »  2NH, 

5) 

AC-i-  3B,  —  AB44B,C 

C043H,  — CH,-H,0 

6) 

2AB,-‘-C,  —  2ACB, 

2PH,4  0,  -  2P0H, 

7) 

Cj-I-AB,  —  ABjC-tBC 

C1,4CH4  -CH,C14BC1 

8) 

2C,  4AB4  —  AB.C,4  2BC 

2C1,4CH4  —  CBjClj-rSHC) 

9) 

3C,- AB4  —  ABC,-  3BC 

3C1,4CB4—  CI1C1,  — 3BC1 

10) 

4C,  4  AB4  —  AC,-4BC 

4C1,4CB,—  CC14  — 4BC1 

A,  B  and  C  label  the  spheres  with  the  refractive  indices  n4  =  1.30,  «„  =  1.50 
and  n,  =  1.40,  and  the  radii  rA  —  l.Ou,,  r„  =  O.BOu,  and  r(.  =  0.75u,,  respec¬ 
tively,  where  «,  is  an  arbitrary  unit.  The  clusters  are  individuated  by  their 
«  chemical  formulae »  in  terms  of  the  «  elements »  A,  B  and  C,  and  their  struc¬ 
tures  are  built  so  as  to  match  those  of  the  chemical  compounds  listed  under 
the  heading  *  EXAMPLE  ».  Since  the  radii  and  the  refractive  indices  of  the 
spheres  are  assumed  not  to  change  due  to  the  « reactions  »,  changes  of  «  valence  » 
or  *  bonds*  are  not  accounted  for.  In  table  I  the  *  reactions*  are  partitioned 
into  tLr«-c  groups:  the  first  group,  *  reactions*  1-44.  where  only  the  spheres 
A  and  B  are  implied ;  the  second  group,  « reactions  »  5  and  6,  where  also  the 
spheres  C  are  implied;  the  third  group,  that  includes  the  family  of  reac¬ 
tions  •  7-410. 

In  table  II  we  list  the  clusters  as  well  as  the  « reactions  »  in  which  they 
are  involved  according  to  table  I.  We  also  report  the  co-ordinates  of  the  centres 
of  the  spheres  to  complete  the  definition  of  the  structures;  in  this  respect  one 
can  observe  that  the  neighbouring  spheres  touch  each  other. 

Figure  1  shows,  for  p  =  1,  with  i)  =  1  w.  I  for  the  «  reactions  »  included 
in  the  first  of  the  above-mentioned  gtoups.  For  completeness  we  report  in 


). 
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fig.  2  the  quantity 

9  =  z  IZT 

Fig.  1.  -  Plot  of  r  r»  i  for  ji  =  1  for  tin-  «  react ions  »  1-4  of  table  J. 

Fig  2  -  Dot  of  p,  u.  I  for  tin-  elu-teir-  implied  in  the  » re.x  t  inns  »  1-4  of  l 1 
Tlie  dashed  curvep  Trier  to  tLr  •  reagents  *  :ui<]  the  solid  curvo  to  tin-  •rt-iiction 
products  *. 

with  the  sanut  value  of  Tt,  r^.  k.  In  the  above  equation  KaZar,  is  the  absorption 
coefficient  of  a  dispersion  containing  only  clusters  of  the  a-th  kind  and  or  numbers 
the  (dusters  that  occur  in  the  «  reactions*  of  the  group  under  consideration 
(in  this  case  the  first).  Furthermore, 

z:o,=  i  o'r", 

to 

Table  II.  -  Co-ordinaUt  oj  the  ctntrr n  of  th<  rpheres  conrtituling  iht  clusters  ivijditd 
in  «  reaction*  »  of  Ui Idt  1.  The  co  ordinates  refer  to  each  sphere  in  the  same  order 
in  which  it  appears  in  the  name  of  the  cluster.  The  « reactions  »  in  which  each  cluster 
is  implied  according  to  table  I  are  also  listed. 


Cluster  Co-ordinates  Reactions 


B,  (—0.5,0.01(0.5.0.0)  1.2,  3.  4.  5 

A,  ( — 1.0,  0,  0)  (1.0. 0. 0>  4 

Ct  (-  0.75,  0.  C)  (0.75.  0.  0)  6.  7.  8.  9.  10 

AB  (—1.0.0,01(0.5.0.0)  2 

AC  (— 1  0.  0,  0)  (0.75.  0.  0)  5 

BC  (-0.5.  0.  0)  (0.75,  0.  0)  7.  8.  9.  10 

B,  (0.  0.577  35,  0»  <0.5. —0.2S8 67.  0j  (-0.5. —0.28S 67.  0)  1 

AB,  (0.  0.  0)  (—1.5.  0.  0)  (1.5.  0,  0)  2 

AB,  (0.  0.  0)  (1.22873. -U.SS036.  0)  (-  1.22S  73. -0.86936.  0)  3 

B,C  (1.02391. -0.71697.  0)  (-  1.023  91. -0.71697.  0)  5 

40,0.0)  '  -  - 

AB,  (0.  0,0)  (0.  1.5.  0)  (1.29901.  —0.75,  0)  (—1.29901.  —0.75.  0)  3 

AB,  (0.  0.  1.38141)  (0.5. —0.8S8  67.  0)  4.6 

(  —  0.5.  —0.28807,  0  (0.  0.57735.  0. 

AB,  (0.  0.  0)  (  — O.S00O2.  -<i.Sll!io2.  V  A'.j(i2)  5.  7.  .if.  I  ) 

(—0.866  02.  0.860082.  —  0.866 02 1 

(0.86602. —0.866  02, —0.86605)  (0.866  02,0.866  02,0.86602) 

ACB,  (0.  0.  1.384  44)  (0,  0.  3.134  44)  (0,  0.577  35,  0)  6 

(0.5,  —  0.288  67.  0)  (-  0.5,  0.288  67.  0) 

AB,C  (0.  0.  1.38444)  (0,  0.57735,  0)  <0.5. —0.28867,  0)  7 

(— 0.5.  —  0.2S8  67.  0)  (0.0.  3.131  44) 

AB, C,  (0. —0.375.  1.45237)  (0.  0.  0)  (0. —0.375.  2  932  37)  8 

(0.75. —  1.0,  0)  (-0.75,  —  1.0,  0) 

ABC,  (0.  1.52069,  0)  (0,  3.02069,  0)  (0.  0.866  02,  0)  » 

(0.75,  -0.43301.  0)  (—0.75,  —  0.433 01.  0) 

AC,  (0.  0.  0)  (-1.01036, —1.01036.  1.01036)  10 

(-1.01036,  1.01036,-1.010  36) 

(1.01036.  -  1.01036.  -  1.01036)  (1.0)036,  1.0)036,  1.01036)  _ 
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whore  o''"1''  is  the  total  cross- section  of  the  i-th  sphere  belonging  to  a  cluster 
of  t hi*  a-th  kind,  calculated  according  to  the  Mie  thoorv  (**).  Therefore  X  Z  0 
is  the  absorption  coefficient  of  the  diversion  of  the  clusters  of  the  a-th  kind 
when  the  multiple-scattering  processes  among  the  spheres  within  each  cluster 
arc  disregarded.  /\  and  giT,  for  the  reactions  of  the  second  and  of  the  third 
group  are  shown  in  the  pairs  of  fig.  3-G,  respectively.  We  notice  that,  since  all 
the  figures  report  the  result  for  rt  ~  1,  this  index  has  been  everywhere  dropped. 

Fig  3.  -  Plot  of  r  rt.  k  for  p  =  1  for  the  •  reactions  •  5-G  of  table  1. 

Fig.  4.  -  Plot  of  <?,  r».  k  for  the  cluhteie  implied  in  the  •  reactions*  5-G  of  table  1. 

Tbt-  dashed  curves  refer  to  the  •  leagentc  »  and  the  solid  curves  to  the  •reaction 
products  ». 

Fig.  5.  -  Plot  of  r  ••*.  1  for  j>  —  1  for  the  «  reactions  •  7-10  of  table  1. 

Fig.  G  -  Plot  of  g a  rt.  k  for  the  clusters  implied  in  the  •  reactions  »  7-10  of  table  1. 

The  dashed  curves  refer  to  the  ♦  reagents »  and  the  solid  curves  to  the  •  reaction 
products  ». 

Figures  2,  4  and  6  are  not  immediately  comparable  with  fig.  1,  3  and  5, 
respectively,  on  account  of  the  different  normalization;  on  the  other  hand,  the 
chosen  normalizations  render  rr,  independent  of  K  and  glT,  independent  of  A'0, 
although  it  should  be  borne  in  mind  that  our  results  are  valid  only  for  low- 
density  dispersions. 

All  our  calculations  were  actually  performed  for  k  ranging  from  0.002ti7> 
to  2.0U”1,  so  that  the  wave-length  of  the  incident  field  was  allowed  to  decrease 
to  about  the  overall  dimension  of  the  clusters.  However,  in  all  figures,  k  ranges 
only  from  O.21171  to  2.0U71,  the  interval  in  which  the  results  show  a  more- 
striking  /.-dependence.  We  remaik,  however,  that  within  the  whole  k-rauge 
the  individual  spheres  have  no  resonance  scattering  as  a  consequence  of  the 
assumed  radii  and  refractive  indices  <**). 

Before  discussing  the  results,  we  premise  a  few  words  about  their  con¬ 
vergence.  In  this  respect  several  authors  deserved  a  particular  attention  to  the 
efTect  of  the  inclusion  of  higher-order  multipolar  contributions  (*•'*’)  so  that 
it  may  be  u.-i-ful  to  point  out  that  our  approach  (>*. -*)  hr  well  as  that  of  ref.  (“) 
is  rather  different.  Actually,  our  formalism,  and  Gorardy  and  Ausloos’  as  well, 
require  that  the  amplitudes  of  the  electric  and  magnetic  multipole  fields  scat* 


(*•>  J- M.  G£kardt  and  M.  Acsloos:  Svrf.  Sei.,  106.  319  (19S1). 

<”)  F.  Clako:  Solid  SlaU  Commvu.,  49,  129  (1984).  / 

(’•)  0.  1.  Si.vdom,  F.  Borgiiese.  P.  Demi  and  G.  Toscano:  Chemical  Systems 
Laboratory  Special  Publication  ARCSL-SP-81009  (September  1981);  O.  1.  Sixboxi. 
F.  Borchese,  P.  Denti,  R.  Saija  and  G.  Toscano:  Aerosol  Sei.  Tech.,  3.  237  (1984). 
(**)  F.  Borchese,  P.  Dextj,  R.  Saija,  G.  Toscano  and  0. 1.  Sindoki:  J.  Opt.  Sot. 
Am.  A,  1.  163  (1964). 
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tercd  by  a  duster  bo  calculated  simultaneously  as  the  solution  of  a  system  of 
linear  equations  whose  ordor  increases  with  tbo  truncation  value  of  L,  say  L1. 
As  a  consequence  any  calculation  with  Z,T  <  Lu  not  only  neglects  tbo  higher- 
ordor  multipoles,  but  also  gives  inaccurate  vain  or  for  tb<-  remaining  multipolar 
amplitudes,  Indeed  tbo  multiple-scattering  processes  among  tbo  spheres  of 
each  cluster  enhance  tbo  strength  of  multipole  fields  of  higher  ordor  (,,-r  ’* 
so  that  tbo  //-expansions  caunot  be  truncated  to  too  low  a  value,  whereas  too 
high  values  imply  time-consuming  computations  and  enormous  memory  re¬ 
quirements.  As  a  matter  of  fact,  wo  found  that  a  quite  reasonable  Lu  =  4 
ensures  converged  results  for  all  the  clusters  considered  here  with  h  ranging 
from  0.002a-1  to  2.0u“\ 

From  a  general  point  of  view,  fig.  2,  4  and  6  show  clearly  that  the  multiple- 
scattering  processes  amoug  the  spheres  of  a  cluster  cannot  be  disregarded  without 
affecting  dramatically  the  correctness  of  the  results.  Therefore,  it  is  a  very 
reasonable  prediction  that  structural  modifications  of  the  clusters  cause  ap¬ 
preciable  variations  of  the  absorption  of  the  dispersion.  As  we  anticipated, 
this  actually  occurs:  the  F-ratios  are  significantly  different  from  unity  for  k 
in  the  range  from  0.002^7*  to  0.2u7\  In  fig.  1,  3  and  5,  however,  we  see  that, 
when  k  increases,  the  effect  of  the  *  reaction  »  becomes  smaller  and  smaller, 
although  for  almost  all  the  *  reactions*  r  is  significantly  different  from  unity 
at  least  up  to  k  ~  1.0u7’.  In  this  respect  we  remark  that,  in  general,  the  more 
the  structure  of  the  clusters  is  modified  the  more  F  remains  different  from 
unity.  Even  the  curves  of  g,  in  fig.  2,  4  and  6  the  less  differ  from  each  other  the 
more  similar  are  the  structures.  As  an  example,  in  fig.  1  the  curves  for  CO, 
and  SO,  are  indistinguishable  from  each  other.  The  component  spheres  are 
identical  for  the  two  clusters  AB,,  indeed,  and  tbe  structural  differences  are 
too  small  to  give  appreciably  different  sj>ectra  also  on  account  of  tbe  random¬ 
ness  of  tbe  orientations.  We  also  notice  tbe  identity  of  tbe  spectra  of  SH, 
in  fig.  2  and  PH,  in  fig.  4.  According  to  table  II,  in  fact,  both  •  compounds  » 
are  AB,  and  have  tbe  same  piramidal  structure.  These  spectra  differ  from  that 
of  SO,  in  fig.  2,  however,  this  latter  « compound  »  is  also  an  AB,  but,  accord¬ 
ing  to  table  II,  Las  a  planar  structure. 

In  conclusion,  the  cluster  model  accounts  for  the  nonsphcrical  shape  of 
the  scatterers  so  effectively  that  it  renders  the  structural  modifications  quite 
detectable.  In  this  sense,  and  on  account  of  tbe  possible  improvements,  tbe 
cluster  model  should  be  suitable  for  approaching  tbe  experimentally  measured 
properties  of  real  aerosols. 
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•  RIASSUKTO 

In  uii  l.ivuin  )>n-c-i-dciite  o  slato  calcolntn  il  ro<flici<,ii1f  di  a>v-orl>iiiicnt(>  jki  ).*>  j>rojia 
paziono  del  c:i)i)]>o  elt-ltroiiiapiMir-n  attravcrso  una  dis]icrf ione  a  Imssa  dmsita  di  dif- 
fusori  non  sfcnri  nindfllizzandoli  conic  •  clnstci  »  di  sfojc.  In  qu<  sto  art  nolo  ci  aininelle 
clic  i  clusU-r  jiossano  camlnaic  stmt  tin  a  in  inanicra  da  simulai*-  d:  lie  reazioni  chimiclic 
e  si  c.iloolano  )<•  ronse-pneM i  vanazioni  del  oooilicionle  di  ahsorliiinento  della  disjaisione 
m  alruni  casi  sipnifiralivi. 


PART  B 


RESONANCE  SCATTERING  FROM  OBJECTS  OF  HIGH  REFRACTIVE  INDEX 


1.  INTRODUCTION 


The  study  of  resonance  scattering  from  dielectric  spheres  is  well 

(1,2) 

known  in  the  literature,  and  its  matematics  is  well  understood  ’  .  However  this 


effect  actually  occur  only  in  spheres  of  very  high  refractive  index,  a  feature 
which  in  practice  is  not  possessed  by  any  particle  of  real  aerosols.  Nevertheless 
the  resonance  scattering  is  able  to  explain  a  number  of  phenomena  so  that  the  as¬ 
sumption  of  high  refractive  index  has  been  accepted  by  several  workers  even  with¬ 
out  a  clear  physical  justification.  We  recall  in  this  respect  that  the  highest 

refractive  index  turns  out  to  be  9  for  water  droplets  in  the  far  infrared. 

(3) 

Recently,  Gerardy  and  Ausloos  pointed  out  that  resonance  scattering 
is  an  actual  occurrence  for  small  metallic  particles  embedded  in  a  dielectric  ma¬ 
trix.  In  fact,  metals  have  a  frequency  dependent  dielectric  constant  with  a  singu 
larity  at  the  plasmon  frequency.  At  this  frequency  therefore  the  scattering  cross 
section  of  small  metallic  particles  presents  a  sharp  peak  due  to  the  high  value  of 
the  refractive  index.  This  phenomenon  is  well  justfied  within  the  framework  of  so 
lid  state  theory. 

The  situation  is  quite  different  for  dielectric  particles  where  the  re 
sonance  of  the  free  electron  modes  do  not  occur.  Therefore  in  order  to  justify  the 


existenceof  high  refractive  indexes  we  took  a  different  approach  based  on  the  ge¬ 
neral-relativistic  treatment  of  the  dielectric  constant  and  magnetic  permeability 

(4) 

of  a  particle  ir  general  motion.  According  to  a  paper  of  Anderson  and  Ryon  ,  in 
deed,  on  account  of  the  covariance  of  Maxwell  equations  under  holonomic  space-time 
transformations,  it  is  possible  to  write  the  dielectric  tensor  of  a  material  me¬ 
dium  in  a  form  that  takes  account  of  the  state  of  motion  of  both  the  medium  and 


the  observer.  Without  going  into  the  calculational  details  here  we  simply  state 
that  the  case  in  which  we  are  interested  is  clearly  that  of  an  observer  at  rest 


while  the  medium,  i.e.  the  scattering  particle,  is  in  general  accelarated  motion 
However,  this  is  not  the  whole  story,  for  to  be  able  to  impose  the  correct  boun¬ 
dary  conditions  across  the  surface  of  the  particle  we  must  also  consider  the  ca¬ 
se  in  which  both  the  medium  and  the  observer  are  in  the  same  state  of  motion. 
Combining  these  two  cases  through  the  appropriate  space-time  transformation,  al¬ 
lows  us  to  conclude  that  for  certain  states  of  motion  of  the  particles,  an  obser 
ver  at  rest  sees  a  scatterer  of  very  high  refractive  index. 

At  this  stage  the  formalism  of  general  relativity  is  no  further  ne¬ 
cessary  for  the  assumption  of  very  high  refractive  index  can  be  inserted  within 
the  customary  theory  of  electromagnetic  scattering.  However,  on  account  of  the 
above  considerations  the  occurrence  of  resonance  scattering  from  dielectric  par¬ 
ticles  gains  a  physical  justification. 


2.  RESONANCE  SCATTERING  FROM  A  SPHERE 


The  theory  of  resonance  scattering  has  been  introduced  within  the  fra 
mewurk  of  the  formalism  we  developed  to  study  the  scattering  properties  of  cluster 
of  spheres.  As  will  be  seen  later  this  work  has  proved  to  be  more  hard  than  expe¬ 
cted,  so  we  will  give  here  only  an  outline  and  present  some  preliminary  results. 
However  it  is  first  useful  to  summarize  the  features  of  the  resonance  scattering 
from  a  single  sphere.  To  this  end  we  recall  that  the  scattering  properties  of  a 


sphere  of  radius  r  and  refractive  index  n 


1 


embedded  in  a  homogeneous  medium  of  re 


fractive  index  n2>  are  determined  by  the  amplitudes  of  the  scattered  multipoles, 

(2) 

aj  and  bg .  As  shown  in  standard  textbook  these  can  be  cast  into  the  form: 


‘  V-O-MicW-  '*■  T*  (N  &0 


(2a) 


k  ^  '*(<&) 

nrw  Cfi)  ~  ^  (P)  “5i 1**0 

where 

a  =  £?T  t*  rriZ/Xa  -  kr,  >  ^2>  = 


Aebeing  the  avelength  in  vacuo  and  n=n^/n2*  The  Riccati-Bessel  functions, if  A  and 

_  (5) 

are  defined  according  to  Kerker 

Now  it  is  evident  that  if  the  denominators  of  eqs(l)  vanish  for  any 

value  of  A,0  and  for  a  given  S. ,  the  corresponding  a^  and  bg  become  infinite  while 

a^,  and  b.  for  £^£'  become  very  small.  This  phenomenon  is  known  as  resonance  scat 

(2) 

tering  of  order  Z  .  Note  that,  as  shown  by  Van  de  Hulst  ,  at  resonance. 

Owing  to  the  presence  of  the  complex  functions  however,  the  denominators  of 
eos.(l)  cannot  vanish  for  any  real  value  of  <£  and  fi  .  Nevertheless,  provided  kr<!C  1 
and  n>>l  there  exist  complex  zeros  with  a  small  imaginary  part.  When  one  puts  the 
real  part  of  these  zeros  into  eqs(l)  one  gets  values  of  a^  and  b^,  that,  although 
very  large  are  not  infinite.  As  a  matter  of  fact  one  has  to  put  krsO.l  or  less 
and  chose  n  of  order  50. 


3.  RESONANCE  SCATTERING  FROM  A  BINARY  CLUSTER 

The  above  considerations  suggest  to  study  the  effect  of  resonance 
scattering  from  a  cluster  of  spheres  in  order  to  get  information  on  the  optical 
behaviour  of  certain  kind  of  nonspherical  impurities  embedded  in  a  homogeneous 
matrix.  We  started  our  study  with  a  binary  cluster  in  order  to  reduce  the  computa¬ 
tional  effort. 

The  first  stage  of  the  research  was  to  find,  for  a  given  kr,  the  re¬ 
fractive  index  which  produces  resonance  in  the  single  spheres  of  cluster.  Actual_ 
ly  the  search  was  made  for  £=1,3,5  and  for  kr=0.1.  The  results  are  summarized  in 


fig.l  where  the  most  striking  feature  is  the  indistingushability  of  the  curves 
for  different  1  even  for  the  value  of  the  main  peak. The  value  of  the  scattering 
cross  section  at  resonance  is  indicated  on  the  figure.  The  n  of  resonance, 
actually  ^ 45  was  then  introduced  into  the  theory  for  the  scattering  from  a 
cluster  for  several  values  of  the  direction  of  incidence.  In  practice  fig. 2  and 
fig. 3  report  only  the  results  for  broadside  (-9’-90°)  and  head-on  (n5>=0°)  inciden 
ce.  The  most  interesting  feature  of  fig. 2  is  the  apparence  of  staellite  peaks  for 
values  of  n  that  are  not  of  resonance,  while  in  fig. 3  some  of  these  peacks  disap 
pear.  This  disappearence  suggest  the  onset  of  some  kind  of  shadow  effect  for  head 
-on  incidence,  but  the  very  origin  of  the  satellite  peaks  is  rather  puzzling,  as 
they  were  quite  unspected  on  the  basis  of  the  general  theory.  We  make  the  hypote- 
sis  that  the  satellite  peaks  originate  from  the  concurrent  effect  of  two  circums 
tances.  First  of  all  the  fact  that,  as  mentioned  above,  a^  b^  ^ ,  at  resonance; 
second,  the  formalism  we  use  to  calculate  the  scattered  field  from  a  cluster,  in 
eludes  the  coupling  of  all  multipoles,  up  to  order  1,  of  the  constituent  spheres. 
Since  even  the  coupling  of  electric  with  magnetic  multipoles  is  taken  into  accou 
nt,  it  may  well  happen  that  the  sum  of  the  above  two  circumstances  produces  peaks 
in  b.  which  though  not  as  high  as  those  of  a  are  shifted  with  respect  to  these 
latter.  In  any  case  this  is  at  present  only  a  reasonable  hypotesis  whose  reliabi_ 
lity  is  to  be  further  investigated. 

4.  CONCLUSIONS  AND  LINES  OF  FUTURE  PROGRESS. 

The  results  illustrated  in  the  fi  ures  are  rather  interesting  but  la 
ck  a  reliable  explanation.  A  careful  analysis  of  the  mathematical  circumstances 
is,  in  fact,  necessary  in  order  to  understand  correctly  the  origin  of  the  satel 
lite  peaks.  This  can  be  done  performing  the  computations  for  highest  values  of  l 


and  for  clusters  of  more  complicated  structures  than  the  simple  binary  cluster 
consided  thus  far.  Furthermore,  the  possible  drop  in  numerical  accuracy  in  the 
vicinity  of  the  peaks  should  also  be  taken  into  account  as  a  consequence  due  to 
the  core  storage  limitations  it  is  necessary  to  modify  the  present  version  of 
our  programs  so  as  to  allow  treating  big  clusters  and/or  large  values  of  Z.  This 
has  already  been  done  during  our  stay  at  the  CRDC,  last  summer.  The  modified  ver 
sion  of  the  program  is  at  present  under  test  and  we  plan  to  use  it  in  the  near 
f  uture . 

The  details  of  the  structure  and  the  use  of  the  modified  program  will 
be  described  by  prof.  Borghese  in  its  next  report  on  Contract  DAJA37-81-C-0895. 
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FIGURE  CAPTIONS 


Figure  1:  Scattering  cross  section  of  a  sphere  for  kr=0.1  as  a  function  of  the 

refractive  index  n.  The  curves  coincide  for  £=1,3,5  and  the  hight  of 

-3 

the  main  peak  is  3.26x10 

Figure  2:  Plot  of  for  a  binary  cluster  as  a  function  of  n  near  the  resonance 
for  broadside  incidence. 

Figure  3:  Plot  of  O  for  a  binary  cluster  as  a  function  of  n  near  the  resonance 
for  head-on  incidence. 
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